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Abstract. For resistor networks on Z'*' with nearest-neighbor conductances, the effective con- 
ductance in a finite set with a given boundary condition is the the minimum of the Dirichlet en- 
ergy over functions with the prescribed boundary values. For shift-ergodic conductances, hnear 
(Dirichlet) boundary conditions and square boxes, the effective conductance scaled by the vol- 
ume of the box converges to a deterministic limit as the box-size tends to infinity. Here we prove 
that, for i.i.d. conductances with a small ellipticity contrast, also a (non-degenerate) central limit 
theorem holds. The proof is based on the corrector method and the Martingale Central Limit The- 
orem; a key integrability condition is furnished by the Meyers estimate. More general domains, 
boundary conditions and arbitrary ellipticity contrasts will be addressed in a subsequent paper. 



1. Introduction and Main Result 

As is well known, most materials, regardless how pure they may seem at the macroscopic level, 
have a rather complicated microscopic structure. It may then come as a surprise that physical 
phenomena such as heat or electric conduction are described so well using differential equations 
with smooth, sometimes even constant, coefficients. An explanation has been offered by homog- 
enization theory: rapid oscillations at the microscopic level average out, or homogenize, at the 
macroscopic scale. However, this does not mean that the microscopic structure is simply washed 
out. Indeed, while it disappears from the structure of the resulting equations, it remains embedded 
in the values of effective material constants, e.g., the coefficients. 

An illustrative example of a homogenization problem is that of effective conductance. We 
will formulate an instance of this problem in the setting of resistor networks. Let denote the 
fif-dimensional hypercubic lattice and suppose that each unordered nearest-neighbor edge {x,y) is 
assigned a value a^y = ayx € (0,°°) — called the conductance of {x,y). For any A C Z^, let B(A) 
be the edges with at least one endpoint in A. Given an /: Z'' — )• M and a finite A C Z^, let 

Ga(/):= £ a,y[f{y)-fix)]\ (1.1) 

where each pair {x,y) is counted only once. This is the electrostatic (Dirichlet) energy for the 
potential / with Dirichlet boundary condition on the boundary vertices of A. 
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Consider now the square box A^, := [0,L)^ nZ'^. A quantity of prime interest for us is the 

effective conductance, 

Cf{t) := inf{eA,(/) : f{x) = ? -x, Vx G ^A^}, (1.2) 

where ? € M'^ and where dA. are those vertices outside A that have an edge into A. By Kirchhoff 's 
and Ohm's laws (see, e.g., Doyle and Snell |7]), Cf^{t) is the total electric current flowing through 
the network when the boundary vertices are kept at voltage t ■ x. 

For homogeneous resistor networks, i.e., when for all (x,j), the infimum ( 11.21 ) is 

achieved by /(x) •=t - x and so C£^(f) = a\t^L'^( \ +o(l)). A question of (reasonably) practical 
interest is then what happens when the conductances axy are no longer constant, but remain close 
to a constant. In particular, we may assume that they are uniformly elliptic, i.e., 

3A G (0,1), V(x,3;) e]B(Z^): X<a^<\. (1.3) 

A, 

A comparison of 2a with these flt^ 's and the homogeneous case shows that Cf^it) is still of the 
order of \t^L'^ . Moreover, thanks to the choice of the linear boundary condition, by subadditivity 
arguments the limit 

ceff(0:=|im i^Cf(0 (1.4) 

exists almost surely for any ergodic distribution of the conductances. The problem left to resolve 
is thus a computation of the limit value. 

Although Ceff (?) can be computed only in a handful of (periodic) cases, it can be characterized 
in large generality: Suppose that a^y = axy((Si) is a sample from a shift-ergodic law P on the 
product space 11 := 0b(z<') [■^j Va] indexed by edges of Z'^. (Technically, axyioi) = ayx{(o) is then 
the coordinate projection on edge (x,j).) As is well known (Papanicolaou and Varadhan lITTl . 
Kozlov 111 Oil and Kiinnermann [11 ] and the book Jikov, Kozlov and Oleinik |9j ), 

ceff(0= inf fif y aQ,,{(0)\t-x + y,g{(0)f\. (1.5) 

Here E is expectation with respect to P, are the unit coordinate vectors in and 

Vxg{(o) := g o Xx{(o) — g{(o) is the gradient of g in direction of x G with denoting the 
shift by x; i.e., the map on Q. such that ay^{Tx(o) := ax+y,x+7{(o) . The expression in (11.51 ) can be 
interpreted as the Dirichlet energy density — with the spatial average naturally replaced by the 
ensemble average. 

Once the (deterministic) leading-order of C£^(f) has been identified, the next natural question 
is that of fluctuations. It is obvious — e.g., by checking the explicitly computable d = \ case — 
that no universal limit law can be expected for general conductance distributions, but progress 
could perhaps be made for the (physically most appealing) case of i.i.d. conductances. However, 
even here establishing just the order of magnitude of the fluctuations turned out to be an arduous 
task. Indeed, more than a decade ago Wehr [21] showed that Var(Cf*) > cU^ for some c > 
but a corresponding upper bound has been furnished only recently by Gloria and Otto lUl. Both 
of these results contain important technical caveats: Wehr requires continuously distributed axy's 
while Gloria and Otto express their results under a "massive" cutoff. 
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Gloria and Otto drew important ideas from an earlier unpublished note by Naddaf and Spen- 
cer ifTSl where (optimal) upper bounds on the variance are derived for certain correlated conduc- 
tance laws. The main tool of [ 15] is the Meyers estimate (cf Meyers fTT]), to be used heavily in the 
present note as well. Other noteworthy earlier derivations of (suboptimal) variance upper bounds 
include an old paper by Yurinskii ||20l and a more recent paper by Benjamini and Rossignol lUl. 
Closely related to these estimates are recent derivations of quantitative central limit theorems for 
random walk among random conductances and approximations of the limiting diffusivity matrix, 
e.g., Caputo and loffe [6], Bourgeat and Piatnitski [i4J, Boivin [3], Mourrat [14], etc. Incidentally, 
the Meyers estimate is also the key tool in 161 . 

The goal of the present note is to prove that, for i.i.d. conductances which are (deterministi- 
cally) not too far from a constant, the asymptotic law of Cf^{t) is in fact Gaussian. Let .yV{iJL,o^) 
denote the normal random variable with mean /i and variance o^. Then we have: 

Theorem 1.1 Suppose the conductances axy are i.i.d. For each d>\, there is A =X{d) G (0, 1) 
such that the following holds: If (11.31 ) is satisfied ¥-a.s. with this X, then for each t (^W^ there is 
G [0,oo) such that 

Cf{t)-^Cf{t) law 2^ 

1]^! -^(0'^' )■ (1-6) 

Whenever the conductance law is non-degenerate we have o} > Ofor all t ^0. 

The proof also immediately yields: 
Corollary 1.2 Under the conditions ofTheorem U.ll 

J-Var(cf (0) o?, (1.7) 

|Az,| L^oo 

where o} is as in (11.61) . 

A few remarks are in order: 



Remarks 1.3 (1) Notice that (11.61) does not give us much information on the "order expansion" 
of Cf^{t). Indeed, we know that 'ECf^{t) is to the leading order equal to Ceff(f)|A/,| but when 
this order is subtracted, the next-order term is (presumably) of boundary size. In J > 3, this is 
still larger than the typical size of the fluctuations. Notwithstanding, what (11.61 ) does tell us is the 
character of the leading order random term. 

(2) There is in fact a formula for of, see Theorem |2.7| below. which also shows that t i-> af 
is a bi-quadratic (and thus smooth) function of t. However, the formula involves complicated 
conditioning and does not seem very useful for practical computations. 

(3) There is no restriction on the single-conductance law other than (11.31 ). In particular, this 
law can have a non-absolutely continuous part including atoms. Certain technical problems do 
arise at this level of generality; see Section 1231 which, we believe, is of independent interest. 



We prove Theorem 11.11 by reducing it to the Martingale Central Limit Theorem. There are 
two main technical ingredients: homogenization theory (which enables a stationary martingale 
approximation of Cf^{t)) and analytical estimates for finite-volume harmonic coordinates (by 
which we control the errors in the martingale approximation). The restrictions to rectangular 
boxes, linear boundary conditions and small ellipticity contrasts permit us to encapsulate the 
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analytical input into a single step, the Meyers estimate, of Proposition 12.41 and Theorem 14.41 
These restrictions can be relaxed but not without additional arguments not all of which have been 
handled satisfactorily at this time. These are deferred to a follow-up paper. 

We remark that two recent preprints have been brought to our attention at the time this work 
was first announced in conference talks. First, Nolen |16| has established a normal approxi- 
mation to the effective conductance defined over a periodic environment, in the limit when the 
period tends to infinity. Second, in a preprint that was posted at the time of writing the present 
note, Rossignol ifTSll formulates and proves a central limit law for the effective resistance for the 
corresponding problem on a torus. Nolen's defines the problem over continuum, albeit with a 
rather strong assumption on an underlying Gaussian i.i.d. structure. Rossignol's setting is based 
on minimizing the electrostatic energy over currents (rather than potentials) subject to a restric- 
tion on the total current flowing around the torus. By a well known reciprocity relation between 
effective conductance and resistance, these papers appear to address similar problems. 

The present paper differs from both Nolen |il6 ] and Rossignol ifTSl primarily in its emphasis on 
fixed (Dirichlet), as opposed to periodic, boundary conditions. Indeed, a majority of our technical 
work is aimed at controlling the resulting boundary effects. Also the way a Gaussian limit law is 
established is quite different: Nolen appeals to Stein's method, Rossignol uses concentration of 
measure while we invoke the Martingale Central Limit Theorem. A deficiency of the present pa- 
per compared to [16 ] and fTSl is the limitation on ellipticity contrast. Nolen overcomes this by an 
appeal to Gloria and Otto [8|, although this ultimately precludes the most interesting conclusion 
in J = 2. Rossignol's approach appears to work seamlessly for all elliptic product laws. 

While the Gloria-Otto method can be adapted to our situation as well, just as for Nolen iT6l it 
fails to deliver the desired conclusion ind = 2. (The issue is that the method yields bounds on the 
moments of the corrector, which diverge ind = 2, while we need only moments of the gradients 
of the corrector.) The moment bounds thus seem to be a separate technical matter, so for our first 
paper we decided to sacrifice on generality of the distribution and derived the CLT only in the 
simplest, albeit still physically interesting, case. 

2. Key ingredients 

Here we discuss the strategy of the proof of Theorem 11.11 and state its principal ingredients in the 
form of suitable propositions. The actual proofs begin in Section |3] 

2.1 Martingale approximation. 

A standard way to control fluctuations of a function of i.i.d. random variables is by way of a 
martingale approximation. Let us order the random variables {axy: {x,y) G B(Ai)} in any (for 
now) convenient way and let to be the a-algebra generated by the first k of them. (Since we 
only aim at a distributional convergence, the a-algebras may depend on L.) Then 



|B{Az,) 



(2.1) 



where 




(2.2) 
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Obviously, the quantity is a martingale increment. In order to show distributional convergence 
to jVi^^o^^, it suffices to verify the (Lindenberg-Feller-type of) conditions of the Martingale 
Central Limit Theorem due to Brown |j5| : 
(1) There exists € [0,°°) such that 

1 |b(Al)I 

£ E(Zf|^,_i) (2.3) 



in probability, and 
(2) for each £ > 0, 



1 



A. 



|B(Az,)| 

I IE(2fl{|z,|>.|A.r/.}|^.-i) -^0 (2.4) 



k=\ 

in probability. 

The sums on the left suggest invoking the Spatial Ergodic Theorem, but for that we would need 
to ensure that the individual terms in the sum are (at least approximated by) functions that are 
stationary with respect to shifts of This necessitates the following additional input: 

(i) a specific choice of the ordering of the edges, and 

(ii) a more explicit representation for Z|. 

We will now discuss various aspects of these in more detail. 
2.2 Stationary edge ordering. 

Recall that ]B(Z^) denotes the set of all (unordered) edges in Z"'. We will order ]B(Z^) as follows: 
Let < denote the lexicographic ordering of the vertices of 11^ . Explicitly, for ;c = , . . . , jc^) and 
y = (ji, . . . we have x <y ii either x = y ox x^y and there exists / G {1, . . . ,J} such that 
Xj = yj for all j < i and Xi < yi. We will write x^yifxy^y and x ^y. 

For the purpose of defining a stationary ordering of the edges, and also easier notation in some 
calculations that are to follow, we now identify B(Z'') with the set of pairs (x, /), where x G 
and / G { 1 , . . . , t/}, so that (x, /) corresponds to the edge between the vertices x and x + e,. We will 
then write 

f either x ^ v 

(x,/) ^ (3;, 7) if <^ . . (2.5) 

^ ^ ' [or x = jandi<j. 

Again, (x, /) -< (yJ) if (x,/) ^ (yJ) but (x,/) / {yJ)- It is easy to check that ^ is a complete 
order on B(Z'^). A key fact about this ordering is its stationarity with respect to shifts: 

Lemma 2.1 If (x, /) ^ {y, j) then also (x + z, /) ■< {y + z^ j) for all z G H' . 

Proof. This is a trivial consequence of the definition. □ 

Now we proceed to identify the sigma algebras {^k} in the martingale representation above. 
Recall that Q. := 0b(Z'')['^i ^/a] denotes the set of conductance configurations satisfying (11.31 ). 
Writing a for elements of Q. we use axy = axy{(o), for {x,y) G B(Z^), to denote the coordinate 
projection corresponding to edge {x,y). Given L > 1, set := |B(Ai)| and let bi, ... ,bN be the 
enumeration of B(Ai) induced by the ordering of edges ^ defined above. Then we set 

^k:=o{o)b: bT^bk), k = \,...,N, (2.6) 
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with 

^o-=(y{(Ob:b^bi). (2.7) 
By definition J^q is independent of the edges in B(Ai) while determines the entire configu- 
ration in B(Ai). Note also that includes information about edges that are not in B(Ai). This 
will be of importance once we replace Zyt by a random variable that depends on all of ft). 

2.3 An explicit form of martingale increment. 

Having addressed the ordering of the edges, and thus the definition of the a-algebras ^j., we now 
proceed to derive a more explicit form of the quantity from ( I2.2I ). Given g) G il, define the 
operator Lq, on (M or M'^-valued) functions on the lattice via 

(U/)(x):= £ a,y{(0)[f{y)-fix)]. (2.8) 

y. {x,y)eBiZ'') 

This is an elliptic finite-difference operator — a random Laplacian — that arises as the generator 
of the random walk among random conductances {axy{co)} (see, e.g., Biskup [2| for a review of 
these connections). The existence/uniqueness for the associated Dirichlet problem implies that 
for any finite A C Z'^ there is a unique *Fa : O. x (AUdA) ^ such that x i-> ^\{(0,x) obeys 

r Lfo»PA(w,x) =0, xGA, 

< (2.9) 
[*PA(ft),x) =x, x£dA. 

It is then easily checked that /(x) := t ■ ^'a{co,x) is the unique minimizer of / i— )• 2a(/) over all 
functions / with the boundary values /(x) = f • x for x G dA. In particular, we have 

Cf{t) = Q^,{t-^'^,) (2.10) 

for all t G M^. The function x H- *Fa(g),x) will sometimes be referred to as a finite-volume 
harmonic coordinate. (The first line in (12.91 ) justifies this term.) 

The minimum value 2a • ^a) is a non-decreasing, continuous and concave function of 
{axy : {x,y) G B(A)}. Thanks to the uniqueness of the solution to ( 12.91 ). 2a '^a) is also contin- 
uously differentiable in axj,'s with 

■^QA{t-'i'A)=[t-'i'A{(0,y)-t-^>A{(0,x)]\ {x,y)£M{A). (2.11) 

This relation is of fundamental importance for what is to come. 

Abusing the notation slightly, let (Oi, . . . ,(0n, with N := |B(A) |, denote the components of the 
configuration co over B(A) labeled in the order induced by ^ defined above. Let 

qi(Oi,...,(ON)-=QA{t-'i'A) (2.12) 

mark explicitly the dependence of the right-hand side on these variables. The product structure 
of the underlying probability measure then allows us to give a more explicit expression for the 
increment Zk = Zk{(Oi, . . . ,(0k): 



J P(dw() . . . F{d(0'fj) [q{(Oi ,...,(0k, (o[^i ,...,(0^) 



- q{(Oi (Ok-i (0^)1 (2.13) 

/r(Ok ^ 
Jco', oCOfc 
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with the inner integral in Riemann sense. A key point is that the last partial derivative is (modulo 
notational changes) given by (12.1 II) . i.e., is the modulus-squared of the gradient of t ■ *Pa over 
the ^-th edge in B(A) integrated over part of the variables. To see that Zyt is a martingale increment 
note that the Riemann integral changes sign when its limits are interchanged. 

2.4 Input from homogenization theory. 

In order to apply the Spatial Ergodic Theorem to the sums on the left of (I2.3j42.4l ). we will 
substitute for Zyt a quantity that is stationary with respect to the shifts of 1/ . This will be achieved 
by replacing the discrete gradient of *Pa — which by ( 12.1 11) enters as the partial derivative of q in 
the formula for — by the gradient of its infinite- volume counterpart, to be denoted by The 
existence and properties of the latter object are standard: 

Proposition 2.2 (Infinite- volume harmonic coordinate) Suppose the law of the conductances is 
ergodic with respect to the shifts ofl/ and assume (11.31 ) for some A G (0, 1). Then there is a 
function Y- D.xZ'' ^W' such that 

(1) (y L(o-harmonic) L(oyA(G),x) = for all x and ¥-a.e. CO. 

(2) (y is shift covariant) For ¥-a.e. CO we have i//^(g),0) := and 

Y{(0,y) - Yi^^x) = Yi^^xO),}! -x), x,yeZ''. (2.14) 

(3) (y is square integrable) 

e( £ ao^{co)\Y{co,x)n (2.15) 

^ j:=ei,...,erf ^ 

(4) (y is approximately linear) The corrector x{(0,x) := i/a(g),jc) — x satisfies 

U„,5MfMP)^„ (2.16) 

|a:|^~ \x\'^ 

Proof. Properties (1-3) are standard and follow directly from the construction of y (which is done, 
essentially, by showing that a minimizing sequence in ( 11.51 ) converges in a suitable L^-sense; see, 
e.g., Biskup fT. Section 3.2] for a recent account of this). As to (4), a moment's thought reveals 
that it suffices to show this for x of the form nci, where n This follows from the Mean 

Ergodic Theorem, similarly as in 12] Lemma 4.8]. □ 

The replacement of (the gradients of) by y necessitates developing means to quantify the 
resulting error. For this we introduce an L''-norm on functions f:Q.x{AUdA)^ W' by the 
usual formula 

/I X'^" 
l|V/||A,p:= I K\f{co,y)-f{co,x)n . (2.17) 

V l^^l (x,v)eB{A) / 

Analogously, we also introduce a norm on functions (p : Q. x Z'' ^ W' by 

\\V(p\\p:=[ £ E\(p{co,x)-(p{co,0)\"y''. (2.18) 

x=ei,...,erf 

Here we introduced the symbol V/ for an M'^-valued function whose i-th component at x is given 
by Vif{x) := /(.x + e,) — f{x) — abusing our earlier use of this notation. It is reasonably well 
known, albeit perhaps not written down explicitly anywhere, that the gradients of *Pa and y 
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close in || • ||A.2-norm (see, however, Proposition 3.1 of Caputo and loffe @ for a torus version of 
this statement). 

Proposition 2.3 Suppose the law P on conductances {axy} is ergodic with respect to shifts of 7/ 
and obeys (11.31 ) for some A G (0, 1). Then 

||V(»Fa,-v^)|| 0. (2.19) 

As we will elaborate on later (see Remark 1431 ). this is exactly what is needed to establish the 
representation (11.51 ) for the limit value Ceff(f) of the sequence L^'^Cf^{t). However, in order to 
validate the conditions (I2.3H2.4I ) of the Martingale Central Limit Theorem, more than just square 
integrability is required. For this we state and prove: 

Proposition 2.4 (Meyers' estimate) Suppose P is ergodic with respect to shifts. For each d>\, 
there is A = A (J) € (0, 1) such that if (11.31 ) holds P-a.5. with this A, then for some p > 4, 

\\Vy\\p < - (2.20) 

and 

sup||V(»Fa,-va)||. <oo. (2.21) 

L>1 ''^ 

Proposition 12.41 is the sole reason for our restriction on ellipticity contrast. We believe that, 
on the basis of the technology put forward in Gloria and Otto lH, no such restriction should 
be needed. To attest this we note that versions of the above bounds actually hold pointwise 
for a.e. ft) € H satisfying ( 11.31 ); i.e., for norms without the expectation E. In addition, from 
m Proposition 2.1] we in fact know (12.201 ) for all /? G (l,oo) when d >3. A torus version of 
Proposition [2]4] appeared in Theorem 4.1 of Caputo and loffe ||6]|. 

2.5 Perturbed corrector and variance formula. 

Unfortunately, a direct attempt at the substitution of (the gradients of) *Pa by Y (12.131 ) reveals 
another technical obstacle: As (12.131 ) relies on the Fundamental Theorem of Calculus, the re- 
placement of *Fa by Y requires the latter function to be defined for ft) that may lie outside of the 
support of P. This is a problem because y is generally determined by conditions (1-4) in Propo- 
sition [2]2] only on a set of full P-measure. Imposing additional assumptions on P — namely, that 
the single-conductance distribution is supported on an interval with a bounded and non- vanishing 
density — would allow us to replace the Lebesgue integral in ( 12.131 ) by an integral with respect 
to ¥{dQ)k) and thus eliminate this problem. Notwithstanding, we can do much better by invoking 
a rank-one perturbation argument which we describe next. 

Fix an index / G {1, . . . ,<i} and recall the notation Vif{x) := f{x + e,) — f{x). For a vertex 
X and a finite set A C Z'' satisfying x G A or x + e, G A, let g^"* (ftJ,x) be defined by 

g'^i!{(0,x)-' :=inf{eA(/):/(x + eO-/(x) = l,/aA = 0}, (2.22) 

where 0^^ := oo. (In Section[5]we will see that q^^ is also a double gradient of the Green function 
for operator Lg,.) Note that (12.131 ) and (12.111 ) ask us to understand how V,*Pa(w,x) changes when 
the coordinate of co over (x,x + e;) is perturbed. Somewhat surprisingly, this change takes a 
multiplicative form: 
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Proposition 2.5 (Rank-one perturbation) Let Kdl/ be finite andx,y GAbe nearest neighbors; 
y = x + tifor some i € {1, . . . ,d}. For any (0, CO' that agree everywhere except at edge b := {x^y), 

V;»Fa(«',x)= [1-(«^-«,)0J)(«',x)]V;»Fa(«,^). (2.23) 
For the prefactor we alternatively get 

l-(w^-Wfo)0j^(«',x) =exp|- / *dWfo0j^(w,;c)|, (2.24) 

where G) coincides with co except at b, where it equals G);,. In particular, 1 — (ft)^ — C0h)Q^f^ {(o' ,x) 
is bounded away from and oo uniformly in (0 ^ Q. and A C 

It is worthy a note that (12.231 ) is a special case of a more general rank-one perturbation formula; 
cf Lemma 15.11 which may be of independent interest. Incidentally, such formulas have proved 
extremely useful in the analysis of random Schrodinger operators. The A t If^-Vamt of the right- 
hand side can now be controlled uniformly in ft) € 

Proposition 2.6 Suppose ( 11.31 ) holds for some Ag(0,1). ThenK^Q^i^{(0,x) is non-decreasing 
and bounded away from zero and infinity uniformly in A C Z'^' and CO (zQ.. In particular, for all 
(O (zO. and all x € Z'^' the limit 

0W(ftJ,x) := lim Q-i\(0,x) (2.25) 

exists and satisfies 

0»(ftj,x)-' =inf{(2z<^(/): + = 1, |supp(/)| <oo}, (2.26) 

where supp(/) := {x G : f{x) ^ 0}. In particular, (ftJ,x) i-)- 0^'' {(0,x) is stationary in the sense 
that 0(')(T,ftJ,x + z) = 0(')(ftj,x) holds for all (O ^Q. and allx,z G Z''. 

Before we wrap up the outline of the proof of Theorem 11.11 let us formulate a representation 
for the limiting variance from Theorem ll.il For x G Z'^ and / G {1, . . . ,d}, let b denote the 
edge corresponding to the pair (x, /) and let 

h{(o,x,i) := / P(dftJ^) r'dcoh [l-{a)h-(Oh)Q'^'Hco,x)]\ (2.27) 



"b 



where &) is the configuration equal to (O except at b, where it equals (bb- Define the matrix 
Z(x, /) := {Zjk{x,i)} j^k=\,...4 by the quadratic form 

(t,Z{x,i)t) :=E(^/i(-,x,/)|V,(f •i//-)(-,x)p o[cOy: b' <{x,i))^, (2.28) 
where (x, /) represents the edge (x,x + e;) and t G M''. Then we have: 

Theorem 2.7 (Limiting variance) Under the assumptions of Theorem \l.l\ the matrix elements 
ofZ{x,i) are square integrable. In particular, of from Theorem \l.l\ is given by 



d . , 

a,2 = £E((f,Z(0,j» ), fGM'^. (2.29) 

i=\ 
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As an inspection of ( I2.28I ) reveals, the limiting variance is thus a bi-quadratic form in t. Al- 
though concisely written, the expression is not very useful from the practical point of view; par- 
ticularly, due to the unwieldy conditioning in (I2.28I ). The representation using the /j-function also 
adds to this; it is no longer obvious, albeit still true, that 

E(^(f,Z(x,/)0 o{cOy: b' -< (x,/))) =0, (2.30) 

i.e., that {t,Z{x,i)t) is a martingale increment. A question of interest is whether an expression 
can be found for that is more amenable to computations. 

Remark 2.8 Since t Cf^{t) is quadratic in t, the above actually implies that, as L — )• oo, the 
joint law of the random variables 

' l^^^li/ ^' --teR'^ (2.31) 
tends to the law of multivariate Gaussian {G, : t € W'} with 

E{G,) = and ^(GfG,) = £ ]E((f,Z(0,;» (5,Z(0,/)5)) , (2.32) 

i=l 

where Z(0, /) is as in (12.281 ). Naturally, f i— G, is a quadratic form as well. 



2.6 Organization. 

The proofs (and the rest of the paper) are organized as follows. In Section [3] we assemble the in- 
gredients — following the steps outlined in the present section — into the proofs of Theorems ll.il 
and l2.7l and Corollary 1 1.21 In Section|4]we then show that the finite-volume harmonic coordinate 
approximates its full lattice counterpart in an L^-sense as stated in Proposition 12.31 and establish 
the Meyers estimate from Proposition 12. 4 1 A key technical tool is the Calderon-Zygmund regular- 
ity theory and a uniform bound on the triple gradient of the Green function of the simple random 
walk in finite boxes. Finally, in Section |5l we prove Propositions 12.51 and 12.61 dealing with the 
harmonic coordinate over environments perturbed at a single edge. 



3. Proof of the CLT 

In this section we verify the conditions (12.3142.41 ) of the Martingale Central Limit Theorem and 
thus prove Theorems 11.11 and 12.71 All derivations are conditional on Propositions 12.3142.61 the 
proofs of which are postponed to later sections. Throughout we will make use of the following 
simple but useful consequence of Holder's inequality: 

Lemma 3.1 For any p' > p>2, a:= and j8 := ^j^, 

Proof. Apply Holder's inequality to the function / := [V(*FAt ~ V^)l- 1^ 

Assume now the setting developed in Section|2l in particular, the ordering of edges and sigma- 
algebras from Section [Z2l and the martingale increment from (12.21) and its representation 
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( I2.13I) from Section [231 In analogy with equation ( I2.27K we also define 



hA{Q),x,i):= / P(dG)^) / dQ)t[\ - {cOt - (Ob)d^^ {a),x)Y , 



(3.1) 



where b := {x,x + ei) and Q) is the configuration equal to co except at b, where it equals (bb- By 
Proposition 12. 51 we may write the martingale increment as 



Zk = E(h/^^{-,Xk,ik)\Vii^{t-'i'Aj{-,Xk) 



(3.2) 



where x^ and 4 are the vertex and the edge direction corresponding to b^, i.e., b^ = {xk,Xk + 
Recall the notation for Z{x, i) from ( 12.281 ) and note that this is well defined and finite P-a.s. thanks 
to the estimates (l2.20j - l2i2TI ) as well as boundedness of h. Note the dependence of on L. 

Proposition 3.2 (Martingale CLT — first condition) Assume that the premises ( and thus con- 
clusions) of Propositions [0]42. 6\ hold. Then Zy_ € L^(P)/or all k and 



1 



£ E(Z,V,_i) £e (?,z(o,/)0 



d 

L 

(=1 



(3.3) 



in ^-probability and L^(F). 



Proof. Fix f e M^. Thanks to Lemma IZT] and Proposition 12.21 2). for each / G {\,...,d}, the 
collection of conditional expectations 



|e(^(?,Z(x,j>)^ o{(0t: b (;c,j))) : xer'l 



(3.4) 



is stationary with respect to the shifts on Z"^ and, by Proposition [23J uniformly bounded in ( 
Labeling the edges in IB (A/,) according to the order ^, the Spatial Ergodic Theorem yields 



1 



£ K{{t,z{xk,ik)tfl^k-i) £e (f,z(o,00 



k=i 



i=l 



with the limit P-a.s. and in L (P). To see how this relates to our claim, abbreviate 

Ak ■=hAi^{-,Xk,ik)\yk{t ■'^Ai,){-,Xk)\ , 



Bk := h{-,Xk,ik)\yk{t ■ Y){-,Xk)\ , 



and denote 



R 



L.k ■- 



^k-x 



(3.5) 

(3.6) 
(3.7) 

(3.8) 



By (|3^ we have Z^ = 'K{Ak\^k), while (12.281 ) reads {t ,Z{xk,ik)t) = E(Bi Hence, as soon 
as we show that 

1 |B(Az,)| 

£ E(|/?L,,|) 0, (3.9) 



A, 



■L\ k=X 



the claim (13.31) will follow. 
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The proof of ( I3.9I ) will proceed by estimating E|/?i^| which will involve applications of the 
Cauchy-Schwarz inequality (in order to separate terms) and Jensen's inequality (in order to elim- 
inate conditional expectations). First we note 

< (E[(A,-B,)2])'^2(E[(A,+B,)2])'^^ (3.10) 

Writing Ak = Bk + (At - Bu) and noting (a + bf < 2a^ + 2b^ tells us 

E [{Ak + Bkf] < 2E [{Ak - Bkf] + 8E (fif) . (3.11) 

Summing over k and applying Cauchy-Schwarz, we find that 

1 

^1 -1 



£ < Ja(2a + 8j8), (3.12) 

where 



1 1 l»V^^Lj 

(x:=j^^ ZE[iAk-Bkf] and /3 := £ E(b2). (3.13) 

By inspection of ( 13.121 ) we now observe that it suffices to show that jS stays bounded while a 
tends to zero in the limit L — )• oo. 

The boundedness of p follows from (12.201) and the fact that h{-,x,i) is bounded; indeed, these 
yield E(|Bi:p) < ||/i||i|fn|VvA||4 uniformly in k and L. Concerning the terms constituting a, 
using {a + b)^ < 2a^ + 2b^ we first separate terms as 

E[{Ak-Bkf] <2E(^\hA,{-,XkA)f\\yi,{t-^Aj{-,Xk)\^-\yk{t-¥){-,Xk))\^^^' 

+ 2E(^\hA,{-,Xk,ik)-h{-,Xk,ik)\^\yk{t-v)i;^k)\'^)- (3-14) 

Since is uniformly bounded, the average over k of the first term is bounded by a constant times 
the product of (||V*Pa^ ||a^,4 + ||Vva||a^,4)^ and i|V(*PA^ - y)\\al4- ^^^^'^ ^^^'^^ ^^""^ 
L — oo by Proposition 12.41 Proposition 12.31 and Lemma [3TT] (with the choices p -.= 4 and p' > 4 
but sufficiently close to 4). 

For the second term in (13.141 ) we pick p > 4 and use Holder's inequality to get 

1 |B{Al)| 2 

£ E(\hAi^{-,Xk,ik)-h{-,Xk,ik)\ |V,-j(f • i/A)(-,x/t) 



\M k% 



< kl'llVl/^H.^ J^\{\hA,{;Xk,ik)-h{;Xk,ik)f')) \ (3.15) 



k=\ 



where q satisfies Vfy = 1. The norm of || Vi/^Ha^.p is again bounded by Proposition 12.41 as long 
as p is sufficiently close to 4; to apply (12.201) . we need to invoke the stationarity of Vi/a to realize 

For the second term in (13.151 ) we first need to show that for each £ > there is A/^ > 1 so that 
for all wen, 

dist^i(2rf)(x,A£) > A/^ |/jAi(ft),x,/) -/j(g),x,/)| < £. (3.16) 
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For this we use that, thanks to ( I2.27I ). dTT] ) and ( fT3] ). 

\hA{(0,x,i) -h{(0,x,i)\ <C J d(5/, |gj('(w,x) -0'^')(w,x)| (3.17) 

for some constant C = C(A) < oo. To estimate the right-hand side, by the monotonicity of A i-^ 
g^j^{a),x) and its stationarity with respect to shifts, we have 

|0!;)(co,x)-g«(w,x)| < |0j^(T^w,O)-gW(T,-w,O)|, w G (3.18) 

as soon as x + A/v C A. Then (13.161 ) follows from (13.171 ) and the fact that the difference on the 
right-hand side of (13.181 ) converges to zero uniformly in ft) € il. 

We now bound the last term in (13.151) as follows. The terms for which is at least steps away 
from Al are bounded by e thanks to (13.171 ): the sum over the remaining terms is of order NL^^^ 
thanks to the uniform boundedness of /ja — h. Hence, in the limit L — > oo, the second term in 
(13.151) is of order e '^i; taking £ J, shows that a tends to zero as L — oo. Invoking (13.121) . this 
finishes the proof of (13.91 ) and the whole claim. □ 



Proposition 3.3 (Martingale CLT — second condition) Assume that the premises (and thus 
conclusions) of Propositions \Z3\\2.6\ hold. Then for each e > 0, 



|b(Al)I 



1 

k=i 



^k-i ] 0, (3.19) 



in W -probability. 



Proof. This could be proved by strengthening a bit the statement of Proposition [3]2] (from squares 
of the Z's to a slightly higher power), but a direct argument is actually easier. 

First we note that it suffices to show convergence in expectation. Let > 4 be such that the 
statements in Proposition l2.4l hold. By Holder's and Chebyshev's inequalities we have 



Since h\^ is bounded, Jensen's inequality yields 



p-4 
2 



l{|z,|>e|A.|V.}) < (^nAyi) ' HW")- (3-20) 



E(|Z,pV2) <CE(^[E(|V4(f-»FA)(-,x,)p|^,)]'^'^ <CE(|V,-,(f-»FA)(-,x,)|^). (3.21) 



It follows that 

1 |B(Az,)| 

-- £ ^\zu\p|^)<c\tnv^>^x,p■ (3-22) 
i^^^i k=\ 

The right-hand side is bounded uniformly in L. Using this in ( 13.201 ). the claim follows. □ 

We can now finish the proof of our main results: 
Proof of Theorems U. 1 \ and \277\ from Propositions 12. 5112. 61 The distributional convergence in (11.61 ) 



is a direct consequence of the Martingale Central Limit Theorem whose conditions (I2.3H2.4I ) are 
established in Propositions 13.21 and 13.31 The limiting variance is given by the right-hand side 
of (13.31) . in agreement with (12.291 ). It remains to prove that a/ > whenever t ^0 and the law P 
is non-degenerate. 
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Suppose on the contrary that o} = 0. Then for each / we would have E{{t,Z{0,i)tf) = and 
thus (f,Z(0,/)0 = P-a.s. Denoting b := (0,e,), dllTtillS] ) imply that, for P-a.e. (Ot, 



CO',) r debt E ( [1 - - «,)gj) («, 0)] I V,(f • V^)(«, 0) | 

■J a',. ^ 



^(o,,-))=0, (3.23) 



where ^(o.;) := O'(g)/,). Let Q.^ C [A,Va] be the set of G)^ where this holds. Then ^{Q.\) = 1 
and, since P is non-degenerate, Q.\ contains at least two points. The expectation in (13.231) is 
independent of (O^; subtracting the expression for two (generic) choices of (Oy in Q.^ then shows 
that the inner integral must vanish for all (Ob,(ol G Q-i. But (12.241 ) tells us that the prefactor in 
square brackets, and thus the conditional expectation, is non-negative (in fact, it is bounded away 
from zero). Hence, this can only happen when 

V,(f-i//-)(-,0) =0, P-a.s. for all/ = 1,...,^^. (3.24) 

But then Cgff (f) = 0, which cannot hold for f 7^ when (11.31 ) is in force. □ 



Proof of Corollar\ \l .2 { from Propositions\2j\^Z6\ Thanks to (llTtilll ) and Proposition [ill Cf{t) 
is a martingale whose increments, are square integrable. Therefore, 



Var(cf(0)= I E(Zf). (3.25) 



k=i 



But the right-hand side is the expectation of the quantity on the left of (13.31 ). Since the convergence 
in (13.31) occurs in (P), the claim follows. □ 



4. The Meyers estimate 

The goal of this section is to give proofs of Propositions 12.31 and 12.41 The former is a simple 
consequence of the Hilbert-space structure underlying the definition of a harmonic coordinate; the 
latter (to which this section owes its name) is a far less immediate consequence of the Calderon- 
Zygmund regularity theory for singular integral operators. 

4.1 bounds and convergence. 

Recall our notation Lo for the operator in (12.81) . We begin by noting an explicit representation of 
the minimum of f (/) as a function of the (Dirichlet) boundary condition: 

Lemma 4.1 Let A C be finite and fix an co £ Q.. Then there is K : dA x dA -> [0,oo), de- 
pending on A and (0, such that for any h that obeys Lc^h^x) = Ofor xG A, 

QA{h) = ^ £ K{x,y)[h{y)-h{x)]\ (4.1) 

^ x.yedA 

Moreover, K(x,y) = K{y,x)for all x,y G dA and 

£ K{x,y) = £ a,, (4.2) 

vedA zGA 

{x,z)eB{A) 

for all X € dA 
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Proof. "Integrating" by parts we obtain 

Q^{h) = -Y^h{y){l^h){y)+ £ a^[h{y)-h{x)]h{y). 

yeA yedA.xeA (4-3) 

'{x,y)et 



Employing the fact that h is -harmonic, the first sum drops out. For the second sum we recall 
that h{x) = T.zedAPA{x,z)h{z), where pa{x,z) is the discrete Poisson kernel which can be defined 
by pa{x,z) '■ = P^{X^^^ = z) for t^j^ denoting the first exit time from A of the random walk in 
conductances ft). Now set 

K{y,z):= £ axyPA(x,z) (4.4) 
xeA 

{x,y)eM{A] 

and note that I,zedAK{y,z) = I.vGA,(;f,>)eB(A) «rv- It follows that 

£ a^[h{y)-h{x)]hiy)= £ K{y,z)[h{y) - h{z)]hiy). (4.5) 

yedA,xeA y,zEdA 
(x,y)eB(A) 

The representation using the random walk and its reversiblity now imply that K is symmetric. 
Symmetrizing the last sum then yields the result. □ 

Remark 4.2 We note that Lemma 14.11 holds even for vector valued functions; just replace 
[h{y) —h{x)]^ by the norm squared of h{y) —h{x). This applies to several derivations that are 
to follow; a point that we will leave without further comment. 



We can now prove Proposition 12.31 dealing with the convergence of V*Fa to Vi/a in |[ • ||a,2- 
norm, as A := A^, fills up all of . 

Proof of Proposition 12.31 Abbreviate h{x) := ^f{(0,x) — *Fa^(g),jc). The bound (11.31 ) implies 

||V(»Fa,-va)||^^2^}-^e( £ a„\h{y)-h{x)A. (4.6) 

^ l^^^l V(.r,y)eB(Az.) / 

Let / : A U 5 A — )■ M'' be the minimizer of 

infj £ \f{y)-f{x)\\f{z)=x{z)fox?Mz&dKi\. (4.7) 

Since h is the minimizer of the corresponding Dirichlet energy with conductances {ary} and 
boundary condition we get using (11.31 ) 



(4.8) 



£ axy\h{y) - h{x)\^ < £ axy\f{y)-f{x)\^ 

{x,y)m{AL) (.t,y)eB(Az,) 

£ \m-m\\ 



(;f,y>GB(Az,) 

Writing the last sum coordinate-wise and applying Lemma I4TT1 we thus get 

£ a^\h{y)-h{xf <^ £ K{x,y)\x{(0,y)-x{(0,x)\\ (4.9) 

(x,y)eB(AL) ^ X.yedAt 

where the kernel K{x,y) pertains to the homogeneous problem, i.e., the simple random walk. 
Note that these bounds hold for all configurations satisfying (11.3b . 
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By shift covariance and sublinearity of the corrector (cf Proposition \2.2[ 2A)). for each e > 
there is A=A{e) such that 

E{\x{;x)-x{;y)f) <A + e\x-y\\ (4.10) 
Using this and (14.91 ) in (14.61 ) yields 

\\'^i'^^L-V)\\l,2^^Tri E K{x,y){A + e\x-y\^). (4.11) 

But T^yedAt^i^^y) — 1 foJ" ^^^^ ^ ^ d-^L while L.Y.veaAL^(-^i3')l-^ ~3'P i^' Lemma RTT] the 
Dirichlet energy of the function x i-^ x for conductances all equal to 1 . Hence, the last sum in 
(14.111 ) is bounded by A \ dA^ \ + e |B (A^) [ . Taking L^oo and £ 1 finishes the proof. □ 

Remark 4.3 As alluded to in the introduction, the -convergence V^Al — > Vi//^ permits us to 
prove the formula (11.51) for Ceff(f)- The argument is similar to (albeit much easier than) what we 
used in the proof of Proposition 13.21 Indeed, we trivially decompose 

Cf{t) = Qa, [t • »Pa,) = Qa, {t-v) + [Qa, {t ■'i'A,)- Qa, (t-v))- (4. 12) 

The stationarity of the gradients of y and the Spatial Ergodic Theorem imply that for any ergodic 
law P on conductances, P-a.s. and in (P), 

r^QA,it-V) L ao4(o)\t-Yico,x)\A. (4.13) 

It follows from the construction of the harmonic coordinate that expression on the right coin- 
cides with the infimum in (11.51 ). (There is no gradient on the right-hand side of (14.131 ) because 
Y{(0,0) := 0.) It remains to control the difference on the extreme right of (14.121) . 

Using the quadratic nature of Qa, the ellipticity assumption (11.31) and Cauchy-Schwarz, 

^QA{t-'i'A)-QA{t-¥)\ 
|A| 

<i|f I|v(»fa-i/.)||^2+xI'I'ii^'^i|2|I^(*^'a-v^)L,2. (4.14) 

By Proposition 12.31 — which holds for any shift-ergodic (elliptic) law on conductances — the 
right-hand side tends to zero as A := A^ increases to Z''. Since we know that |Ai|~^C£^(f) is 
bounded and converges almost surely (e.g., by the Subadditive Ergodic Theorem), it converges 
also in (P). We conclude that the limit value Cgff (f) is given by (11.51 ). 

4.2 The Meyers estimate in finite volume. 

Key to the proof of Proposition 12.41 is the Meyers estimate. The term owes its name to Norman 
G. Meyers |[T3l who discovered a bound on L^-continuity (in the right-hand side) of the solutions 
of Poisson equation with second-order elhptic differential operators in divergence from, provided 
the associated coefficients are close to a constant. The technical ingredient underpinning this 
observation is the Calderon-Zygmund regularity theory for certain singular integral operators 
in M'^. (Incidentally, as noted in |[T3l . Meyers' argument is a generalization of earlier work of 
Boyarskii, cf |[T3l ref. 2 and 3] for systems of first-order PDEs and a version of his result was also 
derived, though not published, by Calderon himself; cf Iil3t page 190]). 
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To ease the notation, in addition to ( I2.18I ). we will use the notation ||/j|p also for the canonical 
norm in ^^(A), 

ll/llp:=(i:|/Wr)''", (4.15) 

throughout the rest of this section. 

Let us review the gist of Meyers' argument for functions on . Our notation is inspired by 
that used in Naddaf and Spencer |[T5l and Gloria and Otto L8J. A general form of the second order 
difference operator L in divergence form is 

L:=V*-A-V, (4.16) 

where A = {A;y(x) : /, 7 = 1 , . . . , (i, x G l/^ are x-dependent matrix coefficients, V/(x) is a vector 
whose /-th component is V,/(x) := /(x + e,) — /(x) and V* is its conjugate acting as V*/(x) := 
fix) — /(x — e,). The above L is explicitly given by 

(L/)(x) = £ (A,-,Xx)[/(x + e,-)-/(x)] -A,-Xx-e,-)[/(x + e, -e,-)-/(x-e,-)]). (4.17) 

Now, if A is close to the identity matrix, it makes sense to write 

L = A + V'^-(A-id)-V, (4.18) 

where we noted that the standard lattice Laplacian A corresponds to V* • id • V. This formula can 
be used as a starting point of perturbative arguments. 

Consider a finite set A C Z'' and let g : A U dK — > M^. Let / be a solution to the Poisson 
equation 

-L/ = V'^-g, in A, (4.19) 

with / := on 5 A. Employing (14. 18I ). we can rewrite this as 

-A/ = V'^-[g+(A-id)-V/]. (4.20) 

The function on the right has vanishing total sum over A and hence it lies in the domain of the 
mverse (A)^' of A with zero boundary conditions. Taking this inverse followed by one more 
gradient, and denoting 

^ := V(-A)X'V*, (4.21) 

this equation translates to 

V/ = JrA-[g+(A-id)V/]. (4.22) 

A first noteworthy point is that this is now an autonomous equation for V/. A second point is that, 
if 11=^ lip is the norm of as a map (on vector valued functions) ^^(A) — )• ^''(A) and ||A — id||oo 
is the least a.s. upper bound on the coefficients of A(x) — id, uniform in x, we get 

||V/||p<||jrA||p||A-id|U||V/||p + ||J^^v||p||g||p. (4.23) 

Assuming || J^||p||A — id||oo < 1 this yields 

"^^"/'^ i ■ (4.24) 

1 - ||J£A||p||A-id||c» 

Furthermore, the condition ||.^||p||A — id||oo < 1 ensures the very existence of a unique solution 
V/ to (14.221) via a contraction argument; (14.241) then implies the continuity of g 1— V/ in £p{A). 
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The aforementioned general facts are relevant for us because Ljq is of the form ( I4.16I ). Indeed, 
set A,y(;c) := 5ijajc.x+Ci and note that (14.171 ) reduces to (12.81) . The finite-volume corrector 

Xa{(0,x) := *Fa(w,x) -X (4.25) 

then solves the Poisson equation 

- ^-coXa = y* -g, where g{x) := (a^^a+gi , • . .,ajc^x+cj)- (4.26) 

Thanks to (II. 31 ). this g is bounded uniformly so, in order to have (14.24b for all finite boxes, our 
main concern is the following claim: 

Theorem 4.4 For each p E (1,°°), the operator bounded in £p{Ai), uniformly in L>\. 

Proof of Proposition \2.4\ from Theorem \4.4\ Let p* > 4. Since (in our setting) ||A — id||oo < 
— 1, we may choose A G (0, 1) close enough to one so that sup^>j ||.;^^ ||A — id||„ < 1. 
From the above derivation it follows 

sup ||V;t;Aj|A.,p* < (4.27) 

L>1 

We claim that this implies 

l|Vz||p<-, P<P*. (4.28) 
Indeed, pick a > and note that, for any £ G (0, a), 

£ l{|V;c(.,A-)|>a} < L ^{\VxAA;x)\>a-e}+ £ l{|VzAj-,.v)-V;t(-,.v)|>e} • (4.29) 
xeAl xeAt xeAi 

Taking expectations and dividing by [A^l, the left hand side becomes P(|V;^(-,0)| > a), while 
the second sum on the right can be bounded by £^^||V;^Al ~^X\\a^2^ which tends to zero as 
L — > oo by Proposition 12.31 Applying Chebyshev's inequality to the first sum on the right and 
taking L^oo followed by £ 1 yields 

P(|V;^(.,0)| > a) < -^sup ||V;^aJ|^; (4.30) 

Multiplying by a''^^ and integrating over a > then proves (14.281 ). 

Returning to the claims in Proposition 12.41 inequality (14.281) is a restatement of (12.201 ). Since 
(I4.27l - |4!28] ) imply the uniform boundedness of ||V(;^Az. — X)\\al,p, for each p < p* , Lemma [3?T] 
then shows || V(;|^Al — X)\\al,p ^ 0, as L ^ oo for all p < p*. This proves (12.211) as well. □ 

4.3 Interpolation. 

In the proof of Theorem 14.41 we will follow the classical argument — by and large due to 
Marcinkiewicz — that is spelled out in Chapter 2 (specifically, proof of Theorem 1 in Sec- 
tion 2.2) of Stein's book [19]. The reasoning requires only straightforward adaptations due to 
discrete setting and finite volume, but we still prefer to give a full argument to keep the present 
paper self-contained. A key idea is the use of interpolation between the strong £^-type estimate 
(Lemma 1431) and the weak ^^-type estimate for J^a^ (Lemma l431 ). Both of these of course need 
to hold uniformly in L > 1. 

Lemma 4.5 For any finite A C U^, the £^{A)-norm of J^a satisfies || J^a||2 < 1. 
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Proof. Let M' be a Hilbert space and T a positive self-adjoint, bounded and inveitible operator. 
Then for all /i G 

{h,T-'h) = sup{2{g,h)-{g,Tg)}. (4.31) 

We will apply this to given by the space (of M- valued functions) i^{A), T := £ — A and 
h := V* •/ for some /: A with zero boundary conditions outside A. Then 

(V*-/,(£-A)-iV*./)= sup {2{gX-f)-e{8,g) + {gA8)} 

geC-iA) 

= sup {2(V^,/)-£(g,g)-(Vg,Vg) -(/,/)} +(/,/) 

(4.32) 

= sup {-{Vg-f,Vg-f)} + {f,f) 

gef{A) 

< (/,/), 

where we used that V* is the adjoint of V in the space of M^-valued functions £^{A) and where 
the various inner products have to be interpreted either for M-valued or R'^-valued functions ac- 
cordingly. Taking £ J, 0, the left-hand side becomes (/, J^a • /)■ The claim follows. □ 

The second ingredient turns out to be technically more involved. 

Lemma 4.6 J^Al of weak-type (1-1), uniformly in L > 1. That is, there exists Ki such that, 
for all L> I, f e t (Al) and a > 0, 

|{zGAi: |^,/(z)I>a}|<fi^. (4.33) 

Deferring the proof of this lemma to the next subsection, we now show how this enters into 
the proof of Theorem 14.41 

Proof of Theorem \4~4\ from Lemma 14761 We follow the proof in Stein |19, Theorem 5, page 21]. 
We begin with the case 1 < p < 2. Let / G ^^(A^) and pick a > 0. Let /i := /l{|/|>a} and 
fi :=/l{|/|<a}. Then 

|{zGAt: \,jeAjiz)\>2a}\ < |{zeAi: >a}| 

+ |{zGAl: |^,/2| >«}|. (4.34) 

Lemmas [4.51 and [43] then yield 

\{z£Al: |^A,/(z)|>a}|<fi^+^2^, (4.35) 
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with ^1 and K2 independent of L. Multiplying by and integrating, we infer 

poo 

\\'^Ajrp=P aP-'\{z&hL: \,J^Aj{z)\>a}\da 

^ pT,£{K,aP-^\fiz)\ l{|j|>„} +K2aP-'\f{z)\^ k\f\<a}) da 

r\f{z)\ ^ /-^ (4.36) 

= pKil,\f{z)\ / aP-^da + pK2Y^\f{z)\^ / aP-'da 

^ Jo , J\f(z)\ 



= ^Ii/(^)r+^Ii/(z)i^ 

proving the assertion in the case \ < p <2. 

For p € (2,00), the fact that is obviously symmetric implies that || J^a||p = where q 

is the index dual to p. Hence sup^>i W'^AlWp < °° for P ^ (1)°°)- D 



4.4 Weak type-(l,l) estimate. 



It remains to prove Lemma 1461 The strategy is to represent the operator using a singular kernel 
that has a "nearly £^-integrable" decay. Let GA{x,y) be the Green function (i.e., inverse) of the 
Laplacian A on A with zero boundary condition on dA. 

Lemma 4.7 The operator J(f\ admits the representation 

[^A'/W] = I I[vPVf GA(x,3;)]/y(3;), (4.37) 

VGA 7=1 

where the superscripts on the Ws indicate which of the two variables the operator is acting on. 
Proof. Since both G\ and / vanish outside A, we have 

[jrA-/(x)] =V,(£Ga(.,j)(V^-/)(j))(x) 

^>'GA ^ 

= £ ({G^{x+^^i,y)-G^{x,y))Y^[fj{y-h)-fj{y)]) 

d (4.38) 

= L L {GK{x+^^uy + ^^j)-Gp,{x,y + ^Cj))fj{y) 

d 

-L L {GA{x + ti,y)-G^{x,y))fj{y). 

j=lyeZ'' 

This is exactly the claimed expression. □ 



Crucial for the proof of the weak-type (l,l)-estimate in Lemma I46l is an integrable decay 
estimate on the gradient of the kernel of the operator J^^: 
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Proposition 4.8 There exists C > independent ofL such that 

|vfV(.iVfGA,(x,,)|<^-^ (4.39) 

for all x^y € and i,j,k G {1, • • • ,d}. 

Although (14.391 ) is certainly not unexpected, and perhaps even well-known, we could not find 
an exact reference and therefore provide an independent proof in Section [43] With this estimate 
at hand, we can now turn to the proof of Lemma l46l 



Proof of Lemma H^ from Proposition 14.81 To ease the notation, we will write A := A^ (note that 
all bounds will be uniform in L) and, resorting to components, write for the scalar-to-scalar 
operator with kernel J(fl^'^\x^y) := V^^W^p GA{x,y) for some fixed i,j £ {l,...,d}. For the 
most part, we adapt the arguments in Stein u9l pages 30-33]. 

Take some function / : A — )• M, extended to vanish outside A, and pick a > 0. Consider a 
partition of into cubes of side 3*^, where r is chosen so large that 3^'^^||/||i < a. Naturally, 
each of the cubes in the partition further divides into 3^^ equal-sized sub-cubes of side 3''^^ which 
subdivide further into sub-cubes of side 3''^^, etc. We will now designate these to be either good 
cubes or bad cubes according to the following recipe. All cubes of side 3'' are ex definitio good. 
With Q being one of these sub-cubes of side 3'^^^ we call Q good if 

^I|/(z)|<«, (4.40) 
1^1 zee 

and bad otherwise. For each good cube, we repeat the process of partitioning it into 3'^' equal-size 
sub-cubes and designating each of them to be either good or bad depending on whether (14.401 ) 
holds or not, respectively. The bad cubes are not subdivided further. 

Iterating this process, we obtain a finite set ^ of bad cubes which covers the (bounded) region 
B := \Jq^^ Q. We define G := Z'' \ B, the good region, and note that 

|/(z)|<«, Z€G, (4.41) 

and 

« < ^ L I ^ 2 € (4.42) 

1^1 zee 

where the last inequality is due to the fact that the parent cube of a bad cube is good. Next we 
define the "good" function 



z)-={ ii: (4-43) 



f{z), z G G 

^Lzeafiz), zeQe 

The "bad" function, defined by b := f — g, then satisfies 

b{z) =0, zeG, 

'£b{z)=0, Q^m. (4-44) 

zee 

Since ■X'Kf = J^Ag + J^Ab, as soon as 



|{z:|^Ag(z)|>«/2}|<^|^ AND \{z:ljeAbiz)\>»/2}\<^^, (4.45) 
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the desired bound ( I4.33I ) will hold. We will now show these bounds in separate arguments. 

Considering g first, we note that 1 1 g 1 1 2 is bounded by a constant times a 1 1 / 1 1 1 . Indeed, for z G B 
let denote the bad cube containing z- Then 

zeZ-i zeB 

<«I|/(^)| + l(^I/(^) 

zeG zeB yeQ, ' (4 45) 

<a||/||i+3'^a£-i-£|/(z)| 

zeB ItJzl 
<(3' + l)«||/||i 

by using (14.41b on G and (14.421 ) on B. By Chebychev's inequality and Lemma 1431 

|{,:|^,,(,)|>„,|<^<Pl±lMm^. ,4.47) 

Note that this yields an estimate that is uniform in A := because || J^||2 < 1 by Lemma 1431 

Let us turn to the estimate in (14.451) concerning h. Let {Qk'- ^ = 1, . . . , |^|} be an enumeration 
of the bad cubes and let := hlq^, be the restriction of h onto Q]^. Abusing the notation to the 
point where we write Jt\{x^y) for the kernel governing ^1^, from ( 14.441 ) we then have 

^A^^(z) = I {^K{z.y)-^K{z.yk)\h(y), (4.48) 

where is the center of (remember that all cubes are odd-sized). Let denote the cube 
centered at j]^ but of three-times the size — i.e., is the union of with the adjacent 3'^ — 1 
cubes of the same side. The bound now proceeds depending on whether z € 2<: or not. 

For z^ Qk, the distance between z and any 3^ G 2*: is proportional to the distance between z 
wdyu- Proposition 14. 8 1 thus implies 

Ka(z,3')-^a(z,w)| <C^^^^^%^, z^^Qk- (4.49) 
Moreover, thanks to (14.431 ). 

I Hy)\ < I {\f{y)\ + \s{y)\) <2 ^ \f{y)\. (4.50) 

yeQk yeQk yeQk 



Using these in (14.481 ) yields 



\^Mz)\<C^^^l^m\. (4.51) 



\^-ykl yea 



Summing over all z^Qk and taking into account that |z — j/tl > diam(2/:) for z G Qk, we conclude 

I W.(z)|<cdiam(G,) I [/(y)! I uTTprr 
<c I 1/(^)1 

yeQk 
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for some constant C. Setting B : = 2/t ^rid summing over k, we obtain 



I m{z)\<cY^\f{y)\<C\\fh, (4.53) 

zga\b ves 

which by an application of Chebychev's inequality yields 

|{ZGA\B: |jrAKz)|>«}|<^^. (4.54) 

i.e., a bound of the desired form. 

To finish the proof, we still need to take care of z G B. Here we get (and this is the only step 
where we are forced to settle on weak-type, estimates). 

It € B; m{z)\ > a}| < |B| < 3''£|eil 
The bound (l433l) then follows by combining (14471 ). (l434l) and (14351) . □ 



4.5 Triple gradient of finite-volume Green's function. 

In order to finish the proof of Theorem 14.41 we still need to establish the decay estimate in Propo- 
sition |48l This will be done by invoking a corresponding bound in the full lattice and reducing it 
onto a box by reflection arguments. (This is the sole reason why we restrict to rectangular boxes; 
more general domains require considerably more sophisticated methods.) 

For £ > 0, let denote the Green function associated with the discrete Laplacian A on 
with killing rate £ > 0, i.e., G^(-, •) is the kernel of the bounded operator (£ — A)^' on 
This function admits the probabilistic representation 

where X is the simple random walk and is the law of X started at x. This function depends 
only on the difference of its arguments, so we will interchangeably write G^{x,y) = G^{x — y). 
We now claim: 

Lemma 4.9 There exists C > such that, for all £ > 0, all i, j ,k ^ {\ , . . . d} and all x^O, 

|V,VyV,G^(^)| <^pj. (4.57) 

Sketch of proof. This is a mere extension (by adding one more gradient) of the estimates from in 
Lawler llT2l Theorem 1.5.5]. (Strictly speaking, this theorem is only for the transient dimensions 
but, thanks to £ > 0, the same proofs would apply here.) The main idea is to use translation 
invariance of the simple random walk to write G^ (x) as a Fourier integral and then control the 
gradients thereof under the integral sign. We leave the details as an exercise to the reader. □ 

We now state and prove a stronger form of Proposition 14. 8 1 
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Lemma 4.10 There exists C > such that, for all L > I, e > and arbitrary /, 7, ^ € { 1 , . . . <i}, 

|vfV;.^VfGi(x,,)|<— ^ (4.58) 

for all x,y ^ A. and all ij,k G {1, . . . ,d}. Here, the superscripts on the operators indicate the 
variable the operator is acting on. 

Proof. Throughout, we fix L G N and denote A := A^. The proof is based on the Reflection 
Principle for the simple random walk on Z'^. To start, denote Aq := Z''^' x N (abusing our earlier 
notation), write X^"^'' for the d-th component of X and let To := inf{^ > 0: xjf^ = 0}. For y ^ Aq 
with components y = {yi, ■ ■ ■ ,yd), put ro{y) := {y —yd)- The Green function on Aq with 
zero boundary condition is given by 

00 

G%,M^y) = Y.{l + e)-'-'P'{Xt=y,To > k). (4.59) 

The Reflection Principle tells us that, for x,y G Aq, 

{Xu = y,ro<k)= {Xk = ro (y) ,ro<k)= {Xt = ro (y)) (4.60) 

and so 

00 

G%^ix,y) = G'ix,y) - £ (1 + £)-^-ip-^(X, = Tq < k) 



k=0 



.G'{x,y)-l^il + e)-'-'P%X, = roiy)) 



(4.61) 



k=0 



= G'{x,y)-G'{xMy))- 
This holds for all x,y ^ Aq and extends even to 3^ G Aq U dAo, as is easy to check. 

Next, consider Ai := Z''^^ x {0, . . . ,L} and set Tl := mf{k > 0: xjf^ = L}. In analogy with 
^o{y) we define ri(y) := (yi , . . . , 2L — y^). The Reflection Principle again yields 

P" {Xk =y,To>k,TL<k)=P' {Xk = rt (y) ,To>k) (4.62) 

and so 

CO 

Gm i^^y) = L (1 + er'^'P'iXk =y,To> k, tl > k) 

k=0 

= Glix,y) -l^{l+e)-'-'F'{X,=y,To >k,XL< k) ^^"^^^ 

k=0 

= G%^{x,y)-Gl{xMy))- 
In conjunction with (I4.61I ). we thus obtain 

Gl{x,y)=Gl{x,y)-Gl{xMy)) 

= G' {x,y) - G' {x, roiy)) - G' (x, ^(y)) + G' (x, ro o rtiy)) 



(4.64) 



for all X G Ai and j G Ai U dAi . 

Proceeding by induction along coordinate directions, we may progressively confine all coordi- 
nates of the random walk and obtain a representation for the corresponding Green function in A 
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in terms of its full-lattice counterpart. Indeed, 

Arf:={0,...,L}^ 



A/ 
Ao 



x{0,...,L}', i=l,...,d-l, (4.65) 



and denoting by , resp., , i = I,... ,d, the reflections in the planes with the i-th coordinate 
equal to zero, resp., L, we get in analogy with ( 14.641 ) the relation 

G%ix,y) =G%Jx,y)-G%Jx,rl;\y)) 

r\ 1^ (4.66) 

Solving the recursion, this leads to 

Gl{x,y)= Y,{-^f^G'{x,Rio...oR,{y)), (4.67) 

where D\ := (8)ti{id 

'''o^'''l^''o^ °f possible reflections that may occur and 

\R\:= Lf^i (Ip _ (,) + 1 (,)) has the same parity as the number of reflections involved. 
From (14.671) and Lemma 14^91 we thus obtain 



|vf)vWvf Gi(x,3.)| < £ \V^VjVlG^{x-R{y))\ < £ (4.68) 

Since, |a; — > \x — y\ when x,y € A, and |fH| = 4^, the result follows. □ 

We are now ready to complete the proof of Theorem I4.4t 

Proof of Proposition 14. 81 Although the £ | limit of G^ exists only in d >3, for gradients we 
have VG{x,y) = limeioVG^{x,y) in all J > 1. Since the bound in Lemma |4 . 1 01 holds uniformly 
in £ > 0, we get the claim in all J > 1. □ 



5. Perturbed harmonic coordinate 
In this section we will prove Propositions I2.5l and l2.6l Abandoning our earlier notation, let 

GA{x,y;(o) = {-La,)-\x,y) (5.1) 

denote the Green function in A with Dirichlet boundary condition for conductance configura- 
tion CO. (Thus, the simple-random walk Green function from Section |4] corresponds to CO := 1.) 
The Green function is the fundamental solution to the Poisson equation, i.e., 

(-loiGA{x,z,co) = 5jc{z.) ifzeA, ^^^^ 

\ GAix,z,co) =0, ifze^A, 



where 5v(z) is the Kronecker delta. Note that Ga is defined for all co & D.. The solution to (15.21 ) 
is naturally symmetric, 

GA{x,y;(o) = GA(y,x;co), x,yeA, (5.3) 

and so we can extend it to a function on A U dA by setting Ga(x, • ; G)) = whenever x G dA. Here 
is a generalized form of the representation (12.23b : 
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Lemma 5.1 (Rank-one perturbation) For a finite K<Zl/ let x,y ^ Kbe nearest neighbors. For 
any (0, co' such that co^ = (Oh except at b := {x,y), and any z e A U dA, 

»Pa(«',z)-»I'a(«,z) 

= -{(0^-co,ry.)[GA{z,y;(o')-GA{z,x;(o')] [»Fa(«,);) -»Fa(w,x)] . (5.4) 

Proof. Suppose ft), ft)' G Q are such that co' equals (O except at the edge b := {x,y), where co^ := 
(Ob + e. Define the function : A U — > M'' by 

<J>a(z) := *Pa(«,z) - £ [GA{z,y; co') - Ga{z,x- co')] ['i'Aico,y) - »I'a(«,^)] • (5.5) 

We claim that 

L^,<I>A = in A. (5.6) 

Since for z G 5A we have <I>a(z) = ^a{(»,z) = z, this will imply <I>a(-) = ^Ai(o', •) thanks to the 
uniqueness of the solution of the Dirichlet problem. 
In order to show (15. 6I ). we first use (15.2145.31 ) to get 

U'<I>a(z) = U'»FA(ft),z)-£[5,(z)-5,(z)] ['i>A{(o,y)-^A{(o,x)]. (5.7) 

To deal with the term L(i)'^a{(0,z), we think of of Lg,/ as a matrix of dimension |A[. For its 
coefficients Lo){z,z') := (5^, Lco5^/)£2(a) we obtain 

lcn'{z,z') = U(z,z') + e [5y.{z) - 5,{z)] [5y{z') - • (5-8) 
Using that Lo*Pa(w,z) = for z G A, we now readily confirm (15.61) . □ 

Proof of Proposition 12. 51 Sety :=.x + e, and denote V,/(z) := /(z + e,) — f{z). Lemma l5?T] shows 

yi^>A{(0\x) = [l - (ftJ^-ftJi)Vj'Vf' GA(x,x,ftj')] V,-»FA(ftJ,x), (5.9) 

where the superindices on V indicate which variable is the operator acting on. To prove the claim 
we need to show 

[vP Vf )Ga(x,x, co)Y'= inf { QA{f) : f{y) - fix) = 1, /^^ = 0} , (5. 10) 

where the conductances in Qa correspond to (O. For this, let / be the minimizer of the right-hand 
side. The method of Largrange multipliers shows 

-U/(z) = a[5y(z)-5,(z)]. (5.11) 

Thanks to ( 15.21 ). this is solved by 

f{z) = a[GA{y.z:,(o)-GA{x,z;co)\ = avS''GA(x,z;ftj) (5.12) 

which in light of the constraint f{y) — f{x) = 1 gives a = [yfhf^GA{x,x,(o)]-\ Since also 
Oa(/) = (/, — Lo/ )£2(A )> (15.111 ) gives 2a(/) = « and so (15.101) holds. The correspondence (12.231) 
then follows from ( 15.9145. lOT) : the identity (12.241 ) results by differentiation of the left-hand side 
with respect to ft)^. □ 

Finally, it remains to establish the limit (12.251) . including all of its stated properties: 



Proof of Proposition \2.6\ Thanks to ellipticity restriction (II. 31 ). we have a bound on this quantity 
in terms of the lattice Laplacian. This shows that, for some c = c(A) G (0, 1), 

c < VpVf'GA(x,x,ftj') < 1/c (5.13) 
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uniformly in A. Moreover, A Vj V- Ga(x,x, co') is obviously non-decreasing in A and so the 
limit exists. The formula (12.26b and the claimed stationarity then follow as well. □ 
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